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NATIOKAT. AmTBQRT-COKaCEmiilE FOE AERONAUTICS 


TECHNICAL NOTE NO. 1429 

lEE EFFECTS OF YAWING THIN POINTED 
WINGS AT SUPERSONIC SPEEDS 
By John C. Evvard 


SUWMAEY 

An approximate relation is derived for the surface velocity 
potential of thin pointed -wings at supersonic speeds -when they are 
contained within the Mach cone from the vertex. This relation is 
applied to obtain the pleasure distributions, the lift and drag 
coefficients, the center of pressure, and the rolling moments as 
a function of angle of yaw for the delta wing. Theoretical cxirvos 
are presented for a Mach ni'mber of to illustrate the relations. 


INTRODUCTION 

The linearized equation of compressible flow at supersonic 
speeds has been applied by Stewart (reference l) and Bro-wn (refer- 
ence 2) to obtain the lift distribution of a thin delta wing with- 
out yaw. Special cases of yawed delta wings have been solved by 
Hayes of North American Aviation, Inc. All these solutions are 
simplified by the fact that the flow is conical. This simplifica- 
tion is inapplicable if the leahing edges are curved. 

A general method has been developed (the basic principle of 
which is presented in reference 3) for ob-balning the lift distri- 
bution of thin wings of arbitrary plan form and profile, Althou^ 
the method may be applied in principle to obtain the effects of 
yawing the thin pointed wing, the practical e-valuatlon of the 
integrals appears to be difficult even by numerical methods. Never- 
theless, an approximate solution to the problem may be obtained 
that will indicate the effects of the plan-form leading-edge curva- 
ture and of yaw on the aerodynamic coefficients of the wing. 

An approximate solution for the surface velocity potential of 
thin pointed wings at supersonic speeds when the leading edges are 
Included within the Mach cone from the vertex was developed during 
March 1947 at the NACA Cleveland laboratory and is presented 
herein. The solution is applied to calculate the pressure 
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distriliution, the center of pressure^ and the wave lift and drag 
coefficients of the delta ving as f mictions of yaw angle. 


ANAT.ySIS 


Tlte basis of the derivation for the surface velocity potential 
is to I'oplace the influence of the pointed wing and the flow field 
between the leading edge and the foremost Mach waves by^ an approxi- 
mately equivalent wing surface without external disturbing flow 
fields when the velocity potential at some local point on the wing 
is calculated. The approximations that are applied will be 
described in detail for the case of the delta wing. 


WlLen the leading edges of the delta wing are swept back along 
the Mach lines (fig. 1(a)), there is no external flow field repre- 
senting Interaction between the top and bottom wing surfaces. The 
correct values of the velocity potential at any point on the top 
wing surface may therefore be obtained through the methods of 
reference 4 by an integration over the peirallelo^am area Sy. 
bounded by the leading edge aiid the forward Mach cone from the 
point (x,y)j that is 


where 

cp 

U 

a 


Cp 


Ua 


dt dm 

/ (x-|) - p (y -t) 


( 1 ) 


-’w 


velocity potential at point (x,y) on top wing 
surface 

free -stream velocity 

angle of attack (negative flow -deflect ion angle on 
top wing surface in y «= constant plane) 


or I, y or Ti Cartesian coordinates 

3 r. ^ --- 


M 


V -1 

free -stream Mach number 


(See references 3 or 4.) 
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If the wing is swept behind the Mach cone fi’om the vea’tex but 
■ is so yawed that one of the leading edges coincides with the Mach 
line (fig. 1(b)), the velocity potential will still be correctly ’ - 
obtained by integrating equation (l) over the parallelograia aa’ea 
defined by the forward Mach lines from the point (x,y) and the 
wing leading edges . 

According to the methods of reference 3, the contribution to 
the potential at point (x,y) from the external flow field Sp is 


= - a 


d| dp 




-/(x-l)' 


A = - cp ^ 

2/„ _2> 


P^(y-n^) 


where cp^ ^ is the portion of the potential contributed by the 
wing region S-^^l. 1!he external field Sj) thus effectively dan- 
ce Is the influence of the wing portion S^^l. 

The case including yaw with both leading edges swept behind 
the Mach cone from the vertex may now be considered. (See 
fig. 1(c).) [fhe potential at point (x,y) is influenced by both 
the wing surfaces S^(o+l+2+3) external fields off the 

wing surfaces Sj,(]_^ 2 + 3 + 4 ) , potential is then 


cp 


IJa 

rt 


®w(0+ 1+2+3) 


' d| dn 

y (x-l)^ - P^(y-T))^ 


U 

TC 






®D(l+2+3+4) 


Adt dp 

y (x-l)^ - P^(y-T])^ 


(2) 


where A represents the slope of the stream sheet in the external 
flow field near the x,y plane measured in = constant planes. 
By the methods of reference 3 



g cL| dfi 


(3) 


and 



A 

XdUn 

[ g d^ dq . 





|%»(1|5|4) 


Adt dri 


1^(143) 


ad^ dn 

-/(x-e)^ - p^(y-T))^ 


Substitution of equations (3) and (4) into equation (2) yields 


iF^ 


( 4 ) 


m 1= 


Ua 


p 


r» 

' 

T 


d^ dq. _ 'Od: 


dUn . u 



, 7(1-02 -p2(7-,)2 ' 

■ J 


7(x-£)2 - p 2(,.;,)2 ” 

%3 

'' u 

£ 


%)(3-h4) 


: Xdj du 

y(i-E)2-p2(y^qT^ 

(5) 


Because X and a have the same sign (reference 3), the second and third Integrals of 
equation (5) tend to counteract each other. ?urthennore^ as 1:he sveepback of the leading edges 
approaches the }fech lines frcrai the vertex, the areas Sy 3 and Sd( 34 . 4 ) approach zero and 

Increases. (!The llmitiDg case is shown in fig. 1(a).) The swcnd and third integrals of 

equation (S) may "Uierefore he neglected, provided that 'the sveepback is not too much greater than 
the Mach angle. Within the validity of this approximation, the potential may he calculated from 
the shaded area (fig. 1 (c)). 
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A comparison of the results thus far obtained shows thet tho 
parallelogram repreaents,tion of the thin delta wing is exact fox* 
the two limiting cases (fig. 1 (a) and fig. 1 (h)), hut is approxi-_~ 
mate for the general case (fig. 1(c)).- The error in the approxi- 
mation is.. likely to he gieatest for wings swept considerably behind 
the Mach cone from the vertex and for wings at zero angle of yaw. 
The paralleiogram approximation may likewise he applied when the 
leading edges of the wing are gently curved on the plan f orm, as in 
figure 2 . 

A set of ohliq.ue coordinates (u,v) having axes parallel to 
the Mach waves simplifies the calculation of the velocity poten- 
tial. In this sytem one of the coordinates of a point is the dis- 
tance measured parallel to the coordinate axis from the point to 
tho other coordinate axis. The transformation epuablons from 
Cartesian to ohliq.ue coordinates are - 


u - (i-Pi) 
V = H (S+Pn) 

f = I (t+u) 


n = i (t-u) 

M 


(6) 


Inasmuch as the elemental area in the (u,v) coordinate sys- 


2(3 ... 

tern is - 5 - du dv, transformation of eq^uatlon (l) by relations ( 6 ) 


M 


yields 


Ucx 


du dv 


S 


’w ,0 


y (u,^-u) (t^-t) 


(■ 7 ) 


where u^ and v^ are the ohliq,uo coordinates of the point (x,y). 
Equation (7) becomes (see fig. 2 ) 



cp = 


Ua 

rtM 
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du 






V 


dv 






. - I 


= ^ ['V - ^aC'^w)] 


( 8 ) 


In the special case of the delta wing (fig. 3), the equations for 
the leading edges are ' 


V = Tj_(u) = 


and 


u » U 2 (v) = 1 c 2 T 

where and k 2 are positive constants, 

In this case, equation (0) hecomes 


c? = ^ y(^w“^l <"w-^^2 -^w) 


^ /ix(l-ki) +py(l+ki)][x(l-k2) - 8y(l+k2)] (9) 


Now the pressure coefficient Cp is 

■ G ^ 

- q " U ^ 


( 10 ) 


\diere dp is the local static pressure minus the free -stream static 
pressure and q is the inooinprossible dynamic pressure (p is 

density). Insertion of equatlcm (9) into equation (10) yields 


'V 


H y(i-ki)(i-t2) 



^-^^1 py 
^ ^ 1-k, X 

" + 

l+k2 gy 
^ 1-kp X 


1+^2 &y K , 

i-t k-< o „ 

1 + *- gy; 

[v 

■ 1-^2 ^ \J 

l-k^ I 


(11) 






f i rpiii : 
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By application of aq.uations (6) and -with the aid of figure 3, 


k = 1 - g tan(9+\lj) 
^ 1 + P tan (9+ ■'!/)' 


ko = 1 - P tan(e-\l/) 
“ 1 + 3 tan(9-\{/)' 


( 12 ) 


\ 

■where \|/ is the angle of yaw and 20 is the vertex angle of the 
delta wing. Substitution of eq,uation (12) into eguation (11) 
yields 


C 


P 


4g / tar'Cg+t^) tan(9-\l/) 

^ tan(0+\J/)] Cl+P tan(e-'.k)] 



tan(0+\j/) 

y/'x 

tan(0-\{^ 


. y/^ 

tan(9-\i/) 

^ tanf0+^)' 


I rJ ' 






( 15)4 


This equation indicates the prossxire distr;Lbution on tluer' thin delta 
wing swept behind the Mach angle at angieb of yaw. Th% equation is 
exact if tan 0 or tan(0+4/) or ±^{6-'^) i^.j3qiml^t£L_l/p - The 
greatest error should ocoiur wiien = J >~ and (t^ 0 « 1 /b^ The 
accuracy of the expression Increases with angle "ST yaw as long as 
the wing lies within the Mach cone frora the vortex. 


The case of '^ = 0 
references 1 and 2. For 
reduces to 


C 


P 


4 

it(l+p ban 0) 


has been solved exactly by the authors of 
\}/ = 0, the approximate equation (13) 


-2a tem 0 



(y/y-)^ 

tan^ 0 



0 


(14) 


The bracketed portion of equation (14) gives the same; type varia- 
tion of the pressure coefficient over the wing surface as was 

obtained in references 1 and 2. The factor 4 ^ how- 

jt(l+p tan 0) 

ever, is somewhat different. The uncorrected values of the factor 
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are compared with the true values of the factor (obtained frcai the 
expressions of reference 1 and the table of complete elliptic 


integrals E of the second kind with modulus v l-p^tan^0, refer 
ence 5) In the following table; 


p tan 0 

Uncorrected 
; value of factor, F’ 
i 4 

True value of 

factor, P 
1^ 

It (l+p tan0) 

! E(-yi-p^tun^e) 

1.00 ■ 

' 

0.6366 

0.6366 

.9 

.6701 

.6697 

.8 

.7074 

.7052 

.7 

.7490 

.7432 

• 6 

.7958 

.7835 

.5 

.8488 

.8257 

.4 

.9095 

.8690 

.3 i 

.9794 

.9121 

.2 

1.0610 

.9520 

.1 

1.1575 

.9842 

.0 

1.2732 

1.0000 


Equation (14) (and, consequently, equation (13)) thus appears 
reasonably accurate, even when p tan d is as low as 0,5. 

In numerical calculations involving equation (13), the 
accuracy of the factor 


F' = 


- jt 


ij [l+p tan ( [l+P tein(0-\|;)] 


(15) 


may be Improved by correction to give the true vaJ.ue at \|/ » 0. 
The original F' and the estimated value of F ore plotted 
against and 0 for p = 1 (that is, M = >/2) in figuro 4. 
In practice, for small angles of yaw, F may be set equal to the 
correct value at \!^ = 0 of 


F = 


1 


E( <yi-p^tan^0) 


~ 0 
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(The ezacb value of the factor F has recently "been derived 
by Max A. Eeaslet^ Harvard Lomax, and Arthur L. Jones of the NACA 
Ames laboratory and is, in the notation of this paper, 

F = i / ■ 2^ 

® V P [tan(e+>l/) + tan.(e-‘4/)] 

where E is the co mplete elliptic Integral of the second kind with 
modulus eq.ua! to V 1-G-^ 

Q = tan(e^-^l/) tan(e-\l/) - tan^(e+^^/)3 [l-g^ tan^(0-^ 

3 Ctan(0+\t;) + tan(0-\{/)] 

For pressure-distribution studies, the pressure coefficient is 
conveniently defined with respect to a set of coordinates fixed on 
the wing. The pressure coefficient In tenns of the wing coordinate 
system is derived in the appendix as eq.uatipn. (A3) and is 


S - 


- oF 


Vr 


tan^ 

j 

(f tan 0 \ 

/ 1 ^ ^1 


Vben0-i-tan\l// 

\tan0+tan'4// xq 

tan2\j/ 


/ tan 0 \ 

/ 1 ) n 


\tan0-tanw 

\tan0-tanW x^_ 


// tan 0 \ 

1 ^ 

^1 

j \tan0 -tanxl// 

\tan0-tan4/(/ 



tan 9 \ 

Vtan0+tanW 


( L h 

\tan0+t^W 


(16) 


The wave lift coefficient Cj^ of the yawed wing may be deter- 
mined by Integrating the pressure coefficient over the wing sur- 
face. This integration is perfoimed in the appendix as equa- 
tion (A5) to give 


Cl = 2^ tane _ 

vl - tan^e tan^ 

A comparison of equations (16) and (17) shows that the lift 
coefficient is related to the pressure coefficient at the center of 
the wing throu^ the equation 
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Cl = -It Cp 


(18) 


where Cp is evaluated at = 0. This relation does not contain 

the coefficient F and hence is exact. Variations of— Mach number, 
angles of-attaok and yaw, and vertex half -angle of the delta wing 
swept behind the Mach angle do not alter the relation, A .simple 
experimental pressure measurement will thus give the lift coeffi- 
cient of the delta wing. By the principles of superposition, the 
pressua'e coefficient of the thin delta wing may be obtained aa one- 
haDf the difference of the pressure coefficients on the top and 
bottom surfaces of a symmetrical -prof ile finite-thickness delta 
wing of the same plan form. The wave drag coefficient (neglecting 
leading-edge suction) is a times the lift coefficient or, with 
the us€) of equation (17), is 


Cj) = 


2ita^ F tan8 
V 1 - tan^0 tan^vj/ 


( 19 ) 


Thus both the wave lift and the wave drag coefficients increase 
with yaw (primarily beca\iae of the variation in F, fig. 3), as 

does I*- It should be noted, however, that a is the flow- 
da 

deflection angle of attack of the yawed wing; the relation between 
a and the angle a', ■vdiloh the delta wing (rotated about its 
base) makes with the plane of zero angle of attack, is 


a = a' cos']/ 


( 20 ) 


Inasmuch as the flow is conical on the yawed delta wing, the 
center of pressure lies along a line parallel to the base at a dis- 

tahce Xp (where Xq is the maximum chord) from the vertex. As 

shown In the appendix (equation (AlO)), the wing coordinates of the 
center of pressure are' given as 




= 


3 

^c 

— ^ tan\!/ 





( 21 ) 


The center of pressure therefore lies on a line parallel to 
the X axle passihig througli the midpoint of the trailing edge. 
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The geometric location of the center of pressure is shown in 
figiire 5. Because the coefficient B cancels in the oalculablbn 
of the center of pressure, the expressions obtained are exact. - 

The rolling moment about the bisector of the rertex angle 
(that is, 7 ^ = 0) may be obtained from equations (21) and (17). 

In terms of the maximum chord, the rolling-moment coefficient C, 
is 


n - n - SirotB’ 

- 7- - ■ 

•*-0 




tan 9 tan\]; 

V 1 - tan^S tan^\]/ 


( 22 ) 


The pitching moment about the axis 
all values of yaw. 


x-j^ = ^ Xj, is zero for 


APPLICiffllOKS OF THIORY 

In order to Illustrate the effects of yaw, the pressure coef- 
ficient on a thin delta wing of vertex half -angle ,tan“^ 0.4 at a 
Mach number ,J~Z is presented in figure 6. In the coiaputatione, 
a corrected value of F (dotted curves of fig.' 4) was used. The 
pressure coefficient seems to remain nearly constant along the 
center line and approaches inf ini ty near either leaSing edge. As 
the angle of yaw is increased, the pressiires on the. least -swept 
side of the wing increase whereas the pressiires , on the most -swept 
side decrease. For the case \j/ = 0, the pressure coefficient , on 
the most -swept leading edge (90° sweepbach) becomes 0. (See 
equation (16).) For greater angles of yaw, the pressure coeffi- 
cient is negative over a portion of the wing surface. (The solu- 
tion obtained when the sweepbacJs: of one of the leading edges is 
greater than 90° does not conform to the Eutta-JouhowsM condi- 
tion for so-called subsonic trailing edges. For this reason, the 
■ 4 ^ = 23.2° ctirve is. dashed.) 

Corresponding to the shift in the pressure distribution (and 
center of pressure), a change in the lift and drag coefficients 
occurs. (See equations (17), (18), and (19).) The variation of 
the lift -curve slopes SCl/ot and are presented in 

figure 7. In terms of the f Low-deflection angle of attach a, 
the lift-curve slope increases with angle of yaw. In terms of the 
geometric angle of attack a' (the wing considered to be rotated 
about its base), the lift-curve slope decreases with angle of yaw. 



12 


MCA rtN No . 1429 



The rolling “moment coefficients of- the wing about the bisector 
of the vertex half. -angle are presented as a function of angle of 
yaw In figure 8. These coefficients were calculated from eq.ua- 
tions (23) and (22) for M = -s/t" and d = tan"^ 0.4. The rolling- 
moment coefficient is seen to be proportional to the angle of 
attack and is nearly proportional to the angle of yaw. A con- 
siderable rolling moment- -develops at large angles of yaw. 

RISULTS OF ANALYSIS 

An analysis of the yawed delta wiug at supersonic speeds con- 
fined wibhin the Mach cone from the vertex gave the following 
results ; 


1. Yawing a delta wing shifts the pressure distribution in 
such a manner as to increase the magnitude of- the pressure coef- 
ficient on the least-swept side of the wing and decrease the pres- 
sure coefficient on the most -swept side of the wing. 

2 . The center of pressure of the yawed delta wing lies al.ong 
a line parallel to the free-stream direction passing through the 
center of the delta base and at a distance two-thirds of the maxi- 
m\jm chord from the vertex. A considerable rolling moment may 
therefore be experienced at large angles of yaw. 

3. The lift -curve slope of the delta wing generally depends 
on the Mach number, the angle of sweepbaok, and the angle of yaw. 
The lift -curve slope of a delta wing of semi-vortex angle 

tan~l 0.4 at a Mach number of r/~2 either increases or decreases 
with angle of yaw, depending upon whether the flow-deflection 
angle of attack or the geometric angle of attack is used. The 
ratio of the two lift-curve slopes is the cosine of the yaw angle. 

4. The absolute value of the ratio of lift- coefficient to 
pressure coefficient at the center— of the yawed thin flat-plate 
delta wing is found to be jt . 


Fli^t Propulsion Beseai:oh Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, July 10, 1947. 
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APPENDIX - CAICUL/VTION OF AIEODYNAMIC COEPFICIEflTS OF . .. 

\ 

YAWED DELTA WING 

A set of Cartesian coordinates (3ci,yj_) defined on the vlng 

surface is convenient for calculating the lift coefficient of the 
thin delta "wing. The transformation eq.uations are 


X = cosv{f + sln^r 


y = y-^ co3\{/ - x-j_ slnii/ 

The auantity y/x then becomes 


y _ 

X “ 


— - tan\|/ 

xi- 


tan'll/ + 1 

^1 


(Al) 


(A2) 


The pressure coefficient^ equation (13), then becomes 




oF 


tan^g - tan^ 
^1 - tan^g tan^\j> 


[/ tan g \ 


■>^1 

Vtang+tanw 

Vtang+tanw 

3:1 

/ tan g N 

( ^ ^ 

5^1 

Vtang -tan\|// 

\tang-tanil// 

^1 


/ tan g > 

_ / 1 

/ Ktang-tan\y) 

\tang-tan\Jf/ ^1 


tan g 
Vtang+tan\l// 


Wang+tan\l/} 


yi 


(A3) 


Because the flow is conical, the lift may be evaluated 

^ dy-, 

triangular infinitesimal area dS = — where Xq 

2 

mum chord. The coordinates and the infinitesimal area 
in the following sketch; 


by using a 
is the maxi- 
are shown 
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„ Zitas' tan© 

~ ■ ■ — -.■...-..I I .. ■ .1 — — 

•J l-tan^© tan^ti/ 


(A5) 


provided that 

js±\J/j < tan~^ i 

The center of presaure 'will lie along the mean 

chord or 5oT = — from the vertex. The other coordinate of the 
a. 3 c 

center of pressure may he obtained from an evaluation of the 
expression. • ' 


ptan© 



-tan© 


The integral in the numerator may he evaluated hy eguations 112, 
113, and 172 of reference 6 to give (for the case \|/ < ©) 


P tan 0 




-tan © 


y tan 0 \ ^ ^1 

\tan0-rtQnw ■\tan©ftan\t/// ^ 

/ tan © \ ^ 1 \ ^1 

\tan©-tan\J// ~ \1:an0-tan\l// 


If tan © \ _ f 1 \ ^1 

+ / \tan©-tan\i// \tan©-tan\j.'y 

\ { tan © \ f 1 \ ^1 

V \tan©+tan\|// Vtan 0 +tan\l/y xi 



V TT 
\^1/ 



It tan^j/ tan^ © 
\/tan^© - tan^\|/ 


(A7) 
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The integral of the denominator may he evaluated hy fomulae 111 
and 113 of reference 6 to' give 


tan 0 


-tan ^ 


Y tan 0 \ 

+ ( ^ A 

^1 

\ tan0+tanil / ) 

V,tan0+tan\{/y 

21 

/ ten 0 ^ 

( ^ -i) 

^1 

\^tan0 -tan'll/ J 

V,tan0-tanw 

^1 


{ tan 0 \ 1 

^ 1 N 

HlH 

!\tan6-tan'ii J ' 

^^tan0-tan\!/y 

If tan 0 \ 1 

^ 1 > 


1 \tan0+tan>lv ' 

^^tan0+tan^l^y 

Xi 


4 2rt tan^ 6 

\ 1/ t/ tan^0 - tan^'}/ 

, (A8) 


Therefore, 


■ (A9) 

5q 


(A simj.lar derivation, with due regard to the sign of the radical, 
shows that situation (A9) also applies when ''1/ > 0 provided that. 
the pressure coefficient may he represented hy equation (A3).) The 
coordinates of the center of pressure are then 


^1 


1^0 





> 


- •=- tanv}/ 
3 J 


(MO) 
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(c) Yawed wing with both leading edges swept behind Mach line. 

Figure I. - Concluded. Para I I e I og ram representation of delta wing. 
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Fig. 8 



Fi gu re 8. ~ Variation of ro I I i ng moment Cj wi th ang 1 e of yaw fo r 

Mach number "'IS*. Angle of attack a. measured in radians. 
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